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Abstract. The formalism of Kohn and Sham uses a specific (model) hamil- 
tonian which highly simplifies the many-electron problem to that of nonintcr- 
acting fermions. The theorem of Hohenberg and Kohn tells us that, for a given 
ground state density, this hamiltonian is unique. In principle, this density can 
be chosen as that of the real, interacting system. To obtain the energy, or 
other properties of the real system, approximations are needed. Working with 
non interacting fermions is an important simplification, but it may be easier to 
produce approximations with different choices of the model hamiltonian. The 
feature that the exact density is (ideally) reproduced can be kept in the newly 
defined fictitious systems. Using model hamiltonians having the same form as 
the physical one, that is, being built of one- and two-body operators, allows 
to approach the physical hamiltonian arbitrarily close, and thus a systematic 
reduction of the approximations. 



1. Introduction 

1.1. General context. The knowledge from first principles of the electronic 
structure of atoms, molecules and solids is contained in the iV-electron Schrodinger 
equation that, within the Born-Oppenhefrner approximation (i.e., at fixed nuclei 
positions R Q ), reads 

(1.1) if(r 1 ,...,r w )*(r 1 oi,...,r i v0-jv) = EV^r^n, r N a N ), 
with 

(1.2) H = f + Wee+Vne, 

N 



(1-3) f = -^XX' 

8=1 

N N 

(1-4) W ee = -^——^-^eiln-r.l) 

N / M \ N 

(1.5) K, - -E Epg )-?>-<*>, 

i—l \a—l 1 1 / 2—1 

where v ne is the external potential due to the M nuclei of charges Z a at positions 
R a , and Hartree atomic units, H = 1 (reduced Planck's constant) m = 1 (electron 
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mass), ao = 1 (Bohr radius), e = 1 (electron charge), have been used. Since elec- 
trons obey the Fermi-Dirac statistics, \& must be antisymmetric under particle ex- 
change, ^(r 1 a- 1 ,...,r i a i ,...,rjaj,...,r N a N ) = -^(ricri, TjCTj, r^i, r N a N ), 
where r denotes the three-dimensional electronic position and a the spin degree of 
freedom (f or J,). In what follows we will be only concerned with the search for the 
ground-state energy, i.e., the lowest eigenvalue Eq of Eq. (|1.1|) . 

The methods that both chemists and physicists have developed to find approx- 
imate solutions of Eq. i|l.l[) can be roughly divided in two large groups: wave- 
function methods (traditional quantum chemistry methods quantum Monte 
Carlo [2]) and density methods (density functional theory |3], and density matrix 
functional theory [4], that is somehow in between the two groups). Simplistically, 
wave-function methods start from an approximation for the antisymmetric, nor- 
malized, iV-electron wave- function, 'J approx, and take advantage of the variational 
principle, 

(1.6) E < So™ = mill (*approx|£|* aP pro X >. 

^approx 

Since 'J approx is a 3N dimensional object, wave-function methods are in general 
computationally expensive when the number of particles increases. 

Density Functional Theory (DFT) uses the electron density n(r) as a basic 
variable, 

(1.7) n(r) = N ^ J |*(rcri, r 2 cr 2 , r N a N )\ 2 dr 2 ...dr N , 

(71...CTJV 

a much simpler quantity to handle, resulting in a low computational cost that allows 
to reach system sizes much larger than those accessible to wave-function methods. 
At given number of electrons TV, the Hohenberg and Kohn theorem [5] tells us 
that the ground-state density n(r) of Eq. (|1.7|) completely determines (except for 
an additive constant) the external potential v ne (r) of Eq. (|1.5(l (for the sake of 
simplicity we only consider physical hamiltonians with a nondegenerate ground 
state). Since the kinetic energy operator of Eq l|1.3l) and the electron-electron 
interaction of Eq. I|1.4f) are the same for all systems, a universal functional F of the 
density n(r) can be defined as |6] 

(1.8) F[n;W ee ,f] = min <*|f + W ee \V), 

where, to keep the connection with what we will do in the next sections, we have 
explictly shown the dependence of F on the electronic interaction W ee and on the 
kinetic energy operator T. The minimum search in Eq. Ijl.8|l is performed over all 
antisymmetric wave-functions \P that yield the density n(r) (by definition n(r) also 
gives, by integration, the number of electrons N). The universal functional F can 
be also defined as a Legendre transform [7] 



(1.9) 



F[n;W ee ,f] = sup jimn(4<|f + iy ee + V\V) - J n(i>(r)drj 



where, as in the rest of this paper, V denotes a local one-body operator of the form 
of Eq. I|1.5fl with v ne (r) replaced by v (r). If the exact form of the functional F was 
known, the variational principle would tell us that 

(1.10) E = min ^F[n-W ee ,f} + J n(r) v ne (r) drj . 
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In practice, because we rely on approximations for F, the energy estimated by 
carrying out the minimization in Eq. (|1.1Q(I can be lower than the exact Eq. 

1.2. Kohn-Sham density functional theory. The Kohn-Sham [g] approach 
to DFT introduces another density functional r s [n], 



(1.11) 



T s [n] ee F[n;0,f] = sup |min($|T + V\$) - J n{v)v{r)dx\ , 



where, as in the rest of this work, $ always denotes the wave-function of a spin- 
i fermionic system with zero electron-electron interaction, i.e., in the majority of 
cases, a single Slater determinant (the fundamental idea of Kohn-Sham has been to 
introduce the fermionic statistic in the construction of T s [n]). Since dealing with 
nonintcracting particles is computationally simple, Kohn and Sham 8 proposed to 
search for that particular noninteracting system which has the same ground-state 
density of the physical one. This defines a model system which is usually called 
Kohn-Sham (KS) system. The difference between F[n;W ee ,T] and T s [n] defines 
the Hartree-exchange-correlation functional -Ehxc [n] , 

(1.12) £ H xcM - F[n;W ee ,f] -T s [n] 

that needs to be approximated. The functional derivative, SEn xc [n]/Sn(r) = 
^Hxc( r )) determines the one-body KS potential «Ks( r ) = w„ e (r)+w H xc( r ) that forces 
the N noninteracting electrons to have the same density of the physical system. The 
KS-DFT relies thus on the assumption that, given a physical ground-state density 
n(r), it is possible to find a noninteracting system which has the same ground-state 
density. 

From the functional BhxcM the classical electrostatic Hartree term is usually 
extracted, 



(1.13) E H [n] = ^J drj dr 



, n(r)n(r') 



|r - r'| 

and the remaining, unknown, part is called exchange-correlation energy, 
(1.14) E xc [n] = F[n; W ee ,f] - T s [n] - E H [n}. 

The Hartree functional En [n] describes the electrons as if they were classical charge 
distributions. It is a simple functional of the density, and yields an important part 
of the energy of a many-electron system, but is nonzero also for a one-electron 
density. This "self-interaction part" of En [n] is cancelled by the exact E xc [n] , but 
this does not occur for most of the current approximate functionals, which suffer 
of the so called "self-interaction error" . 

The success of KS-DFT is mostly due to the fact that simple approximations 
(local-density approximation and generalized gradient corrections) for B X cM an( i 
its functional derivative provide practical estimates of thermodynamical, structural 
and spectroscopic properties of atoms, molecules and solids. However, with the 
current approximations, KS-DFT is still lacking in several aspects, in particular 
it fails to handle near-degeneracy correlation effects (rearrangement of electrons 
within partially filled shells) and to recover long-range van der Waals interaction 
energies. The inaccuracy of KS-DFT stems from our lack of knowledge of E xc [n], 
and much effort is put nowadays in finding new approximations to this term |2] . A 
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trend in the current research is to construct implicit functionals of the density: in 
particular, the exchange-correlation functional is divided into exact exchange 

(1.15) E x [n] = mw ee \$) - E H [n], 

and the remaining correlation energy E c [n] = E xc [n] — E x [n]. The exact exchange 
cancels the self-interaction error of _Eh[^], and is an implicit functional of n(r) 
through the Slater determinant The corresponding Kohn-Sham potential must 
be determined via the optimized effective potential method (OEP) |10| . 

(1.16) Eo = te{($v\f + W ee +V ne \$v)+Ec[n* 1l ]}, 

where is the ground state of the noninteracting hamiltonian, T + V. The con- 
struction of a correlation energy functional E c [n] to be used with exact exchange 
in Eq. (|1.16Jl is still an open problem. 

In this work, we review some basic ideas, results, and open questions of a 
different approach: instead of trying to approximate the KS E xc [n] , we change the 
model system defined by Eq. (|1.11|) . 

1.3. Adiabatic connection formula. Before discussing the choice of differ- 
ent model hamiltonians, we report some equations that will be used in the next 
sections. An exact formula for the functional -EhxcM can be obtained via the adia- 
batic connection formalism |11U12| : by varying a real parameter A, the interaction 
u> A (ri2) between the electrons (we have defined r\ 2 — \v\ — r2 1 to denote pair- 
wise interactions that only depend on the electron-electron distance) is switched 
on continuously from zero to l/?"i2, while the density is kept fixed by an external 
one-body potential V x . Each hamiltoninan H x along this adiabatic connection has 
a ground-state wavefunction that yields, by construction, the same density n(r) 
for each A. If w A=0 = and w x=a = l/ri2, the KS Hartree-exchange-correlation 
energy is given by jlll 112) 

dw x (r 12 ) 
dX ' 

where the spherically and system-averaged pair density (APD) / A (f 12 ) is obtained 
by integrating |^I/ A | 2 over all variables but r 12 = |r 2 — r^|, 

(1-18) / A (r 12 ) = N{N ~ 1] /l* A (ri2,R,r 3 ,...,r Ar )| 2 ^dRdr 3 ...dr w . 

a\ ...ajv 

where R = (r t + r 2 )/2. 

2. Changing the model hamiltonian 

2.1. General considerations. Equation (|l.lll) shows that the KS approach 
to DFT introduces a model hamiltonian in which W ee is set equal to zero and the 
one-body potential, «ks ( r ) i is different from v ne (r) (and is obtained by imposing the 
condition that the ground-state electron density of the model hamiltonian be the 
same of the physical system). Solving the noninteracting KS hamiltonian instead of 
the fully interacting hamiltonian of Eq. (|l.ll) is obviously very practical. The exact 
ground-state energy of the physical system can, in principle, be obtained from the 
KS Slater determinant <& via the functional i5 xc [n]. However, there are cases in 
which the restriction that the model system be noninteracting makes the search for 
an approximate -E xc [n] seem like a daunting task. Consider the simple example of 



ShxcM = I d\ f 
Jo Jo 



(1.17) £hxcM= / d\ dr 12 4irr 2 12 f x (r 12 ) 
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Figure 1. The spherically and system-averaged pair density 
f X (r i 2 ) of Eq. (ji"T%|) for the H 2 molecule in the KS system (A = 0) 
and in the physical system (Coulombic interaction). The inter- 
nuclear equilibrium distance R = 1.4 and the extreme stretched 
molecule at R = 20 are considered. 



the H 2 molecule (two electrons in the field of two nuclei of unitary charge spaced by a 
distance R). In Fig. ^ we report the quantity 47rr^ 2 / A (r 12 ) that enters in Eq. (|1.17|) 
for the KS system and for the physical system (A = and A = a, respectively) 
at two internuclear distances R. As shown by Eq. (|1.17(l . the change in the APD 
f x (r 12) when we switch from the KS system to the physical (fully interacting) 
system determines the functional En xc [n]. If this change is not drastic, we can 
expect that universal approximations for E-n xc [n] can work relatively well. This is, 
e.g., the case of the H 2 molecule at the equilibrium distance R — 1.4 a.u., reported 
in Fig. ^ But when we stretch the molecule (e.g., in the extreme case R = 20 
considered in the same figure), we see that the two APD are completely different, 
which means that the term E-n xc [n\ becomes very important: it has to correct the 
very different nature of the KS wavefunction with respect to the physical one. This 
effect is completely system-dependent and it is thus very difficult to include in a 
universal functional of the density. 

Cases like this appear when we have near-degenerate levels in the Kohn-Sham 
system (in the case of the H 2 molecule described above, the energies of the two KS 
states o g and a u become closer and closer as R increases) . The basic idea reviewed 
in this paper is to remove the constraint that the model system be noninteracting 
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in order to keep a reasonable resemblance between the model wavefuntion and the 
real one. 

2.2. Using a modified electron-electron interaction. Obviously, the model 
system must still be less expensive to solve than the physical one. A possible ap- 
proach is to define a density functional for a "partial" electron-electron interaction 
w(ri2) > 0, smaller than the Coulomb interaction VV12, w(r\2) < l/r±2, but differ- 
ent from zero (we restrict our choice to pairwise interactions that only depend on 
the electron-electron distance). Following the idea of the adiabatic connection of 
Sec. II .HI we can make our partial interaction depend on a real parameter /1 in such 
a way that w^^r 12) — > when fi — ► and w^lr 12) — * l/n.2 when \i tends to some 
positive value a, and define 

(2.1) F»[n] = F [n;W^,T] = sup jmLa<*|T + W fl + V\V) -J n{r)v{r)dr 

With this definition, we see that the density functional f[n] switches from the 
Kohn-Sham one of Eq. to the physical one of Eq. I|1.9ll . 

(2.2) F^°[n] = F[n;0,f]=T s [n], 

(2.3) F^ a [n] = F[n;W ee ,f\. 

In analogy with the KS approach, we can ask that the model system with inter- 
action w^(r 12) have the same density of the physical one. This fixes the external 
one-body potential (that we call v^ir)) in our model hamiltonian . The ground- 
state wavefunction of our model system is denoted it is a multideterminantal 
wavefunction that must be computed with one of the standard methods of quan- 
tum chemistry (configuration interaction, coupled cluster, multi-configuration self- 
consitent field,...). In general, if (i is not too large (i.e., if w M is still much smaller 
than the full Coulomb interaction), few determinants describe ^ M quite accurately, 
so that the computational cost can be kept low. The difference F [n] between the 
physical functional F[n; W ee , T] and the partially-interacting 

(2.4) F>] =F[n;W ee ,f]-F[n;W^,f], 

is what we need to approximate, together with its functional derivative that deter- 
mines v ne (r) — u M (r) . 

How to choose lo^rj^)? Two points are important to determine the partial 
interaction: (i) that we can design reasonable approximations for the correspond- 
ing F M [n] , and (ii) that we can solve adequately the hamiltonian of the model 
system. A convenient choice seems to be a long-ranged interaction, i.e., a w M (ri2) 
that behaves as I/Y12 for large ri2, but that is softer than l/r^ for small ri2- 
The reasons for this choice are (i) short-range correlation effects seem to be more 
transferable from one system to another |13j . and they should thus be more easily 
described in terms of an approximate universal density functional; (ii) the tradi- 
tional wavefunction methods of quantum chemistry can reasonably describe ^ at a 
lower cost than the fully interacting system, because of the smaller interaction and 
because of the absence of the electron-electron cusp. For practical reasons (analytic 
matrix elements for both Gaussians and plane-waves, i.e. the most commonly used 
basis sets in quantum chemistry and solid-state physics, respectively), a common 
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choice is 

(2.5) W ^r 12 ) = ^fciHl, 

where erf(x) is the error function (and we thus have a — oo). Other possibilities for 
w^(r 12) have been also explored: they are all, like Eq. (|2.5|l . arbitrary, and an open 
question is wether there is a way to determine w^ L (r 12) according to some optimal 
criteria. The approach resulting from the choice of Eq. (12.511 is the one that we 
describe more in details in the next Sec. Before doing that, for completness we 
also briefly introduce a different choice of the model hamiltonian. 

2.3. Adding a nonlocal one-body operator. Another way to modify the 
model system is by adding a nonlocal one-body operator, Onl, multiplied by a real, 
positive constant g [141 115] , 

(2.6) F9[n] = F[n;f + g6 NL ,W ee }. 

By setting Onl equal to the projector (here reported for closed-shell systems) onto 
the virtual Kohn-Sham orbitals (f>i (r) , 

N/2 

(2.7) NL (r, r') = <5(r - r') - £ &(r)# (r'), 

i=i 

the functional F 9 [n] of Eq. ((2.6(1 switches from the physical one at g = to the 
Kohn-Sham one of Eq. in the limit g — > 00 (since (^IOnlI^) is always > 0, 

as g — > 00 the model system minimizes its energy by making (^Onl^) = 0, that 
is, by occupying only the KS orbitals). This choice of the model system is not 
further detailed here; the interested reader may find additional information in the 
literature (141 1161 Il5] . 

3. Multideterminantal DFT from a long-range-only interaction 

3.1. The functionals. Following the same steps of Sec. 11.31 we can write an 
exact formula for the functional F^[n] of Eq. (|2.4() in terms of the APD f^(r 12). 
It is convenient to use for the adiabatic connection the same partial interaction 
w^(r '12) of Eq. 1)2.5(1 chosen to determine the model system, 

(3.1) F>]= dp,' dr 12 ^r 2 12 r'(r 12 ) dW ™ 



n JO 



dfi 1 

This formula is identical to Eq. ((2.4(1 , except from the fact that the integration over 
the coupling constant starts from the positive value \i instead than zero. 

As in KS-DFT, the functional F [n] can be divided into an Hartree term, 



r — r 



(3.2) E^[n] = - J dr J dr'n(r)n(r') 

and an exchange-correlation term 

(3.3) XM = F>}-E^[n] 



^(|r-r'|) 



that needs to be approximated. The functional E^. c [n] can, in turn, be divided into 
exchange and correlation in two different ways. We can, in fact, define an exchange 
functional by using the Kohn-Sham determinant 



(3.4) ££ [n] = <*l Wee -W |*>-££ 



n 
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and then define the usual correlation energy functional E c [n] , 

(3.5) E'M^XM-XM, 

but we can also define a multidetcrminantal (md) exchange functional |17| by using 
the wavefunction 

(3.6) X,mM = {^\W ee - W»\*») -E&[n], 
and then a corresponding correlation energy, 

(3-7) Kmd H = Kc N - £ x,md M • 

These two ways of splitting exchange and correlation play a role only if we imple- 
ment an optimized effective potential-like scheme, in analogy with what is usually 
done for the exact- exchange KS-DFT [see Eq. (|l.ltj|) ]. In the case of the mul- 
tideterminantal exchange of Eq. I|3.6[) the problem can be reformulated as |17| 

(3.8) E = mi\(^\f + W ee + V n e\^) +E^ md [n^ l) , 

where is obtained by solving the Schrodinger equation corresponding to the 
hamiltonian H ^ = T + + . This equation is the generalization of Eq. (|1.16(l 
to the case of the multideterminantal model system 



3.2. Exact properties of the functionals. From Eq. i|3.1[l it is possible to 
derive exact properties of the functionals of Eqs. (|3.4|) - (|3.7|l in the fi — ► limit |18L 
[19] and in the fi -> oo limit [181120] . 

In the first case, all the functionals E x) E^ c [Eqs. (|3.4|) and l|3.5[) ] and E xmd> 
E^ md [Eqs. (EH and 1(23] tend to the KS functionals of Sec. O This limit 
can be studied with perturbation theory: one finds that the way in which the 
functionals approach the KS ones depends on wether the system is confined |18j 
(atoms, molecules) or extended |19) . 

The large- /i limit is the most interesting, since, as shown by Eq. (|3.1|l . it always 
lies in the range of /i-values for which we want to construct approximations. In 
this limit, all the functionals of Eqs. (|3.4() - (|3.7|l vanish and, since dw^fr^/d/j, = 
2 e -p r 12 ^ n^.ij) shows that their large-/i behavior is determined by the short- 
range part (small ri 2 ) of the spherically and system-averaged pair density / M (ri 2 ). 
It is possible to show |20| that, when [i — > oo, the small-ri 2 part of / M (ri 2 ) is 
dominated by 



(3-9) T(r 12 ) = /(0) 



1 + 2r 12 pi(fir 12 ) + 



where /(0) is the "on-top value" (zero electron-electron distance) of the physical 
(i.e., corresponding to the Coulomb interaction) /(ri 2 ), and the function pi(y) is 
equal to [2U] 
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Inserting Eq. 1)3.9(1 (and, for the case of Eg[n], the KS j ' A ' _0 (r 12), corresponding to 
the KS determinant <!>) into Eq. l|3.1|l we find, for unpolarized systems |21L 1181 120| 



/ . / s2 ~ I 1 



(3.11) E: W=~ / drn(r)^+0(- 



(3.12) BT°°[n] 



(3.13) 2C~[n] 



(3.14) <^[n] 



/(o)- 


4 J 


/ drn(r) 2 


/(o)- 


2 J 


f drn(r) 2 



-/<o>^>+o 



/(0) 



4V2tt 
3 a* 3 



These equations tell us that (i) if we use the definition of exchange with the KS 
determinant, the functional i^M i s ' f° r l ar S e Mi exactly described by a local func- 
tional of the density; (ii) all the other three functional involve the physical "on-top" 
/(0), which is usually not available (its knowledge would require a very accurate 
calculation for the physical system!). However, it is possible to construct approx- 
imations for /(0). For instance, for many systems the estimate of /(0) from the 
local density approximation (i.e., by transfer from the uniform electron gas model) 
is rather good |13l . Another possibility is to use Eq. I|3.9|) to estimate /(0) from 
/ M (0). This estimate has the advantage of being without self-interaction error 20 . 



3.3. Building approximate functionals. Following the same historical path 
of KS-DFT, the simplest approximation one can think of is the local density (LDA), 



(3.15) < c [n] = J dvn{v) [e xc (n(r)) - e^(»(r))], 

where e xc (n) and £x C ( n ) are the exchange-correlation energy per particle of an elec- 
tron gas of uniform density n with interaction l/ri2 |22| and w^(r 12) |23L I19| . 
respectively. 

Generalized-gradient approximations (GGA) for E^ c have been also constructed 
[24l l25l l26l [27] along similar lines of KS-DFT [28j. In this context, we only men- 
tion a conceptually simple approximation based on the exact properties of the 
previous Sec. 13.21 It consists 24] in a rational interpolation between a given GGA 
KS functional at [i = 0, and zero at fi — > 00. The rational interpolation is con- 
strained to recover the exact large- n behavior of Eqs. (|3.11|I - H3.12|I . using the LDA 
approximation for the physical on-top /(0). 

All these ways of constructing approximate functionals are based on the as- 
sumption of transferability of short-range exchange and correlation effects from 
one system (the uniform electron gas) to the others. Another strategy that came 
out recently, and that seems very well suited to describe short-range correlation 
effects, consists in generating realistic APD / M {r-12) along the adiabatic connection 
[with p! > /j,, to be inserted in Eq. 1)3. ljl ] by solving simple "radial" (unidimensional) 
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equations |29l I3U| for a set of "effective geminals" xj^ (^12) 
1 d 2 £(l + 1 



2 r 12 + —2 +<ff( r 12) 



7*12 dr{ 2 

(3.16) ^^|^f'(n2)| 2 =.r'(r 12 



W (r 12 ) = < V>f (r 12 ) 



The effective interaction u>g ff (ri2) is (in principle) the Lagrange multiplier for the 
exact / M {r\2) (i.e., corresponding to the wavefunction ^ along the adiabatic 
connection; in practice w^ s (r 12) is approximated). In these equations the rule 
for the occupancy #j of the effective geminals is chosen, for spin compensated 
systems, to be Slater-determinant-like: occupancy 1 for even £ (singlet symmetry), 
occupancy 3 for odd I (triplet symmetry), up to N(N — l)/2 occupied geminals. 
This rule has been applied to solve the effective equations l|3.16ll in the uniform 
electron gas at fj, = 00 (Coulombic interaction), with rather accurate results |31U32) 
when combined with simple approximations for the effective interaction potential 
w cS {r 12 ). 

3.4. The calculation of VP M : an illustrative example. As said in the 
previous paragraphs, the model wavefunction is, in most cases, computed with 
one of the traditional wavefunction methods of quantum chemistry. We report here 
an illustrative example of such calculation, and we then give in the next Sec. 01 
an overview of recent results obtained using different methods for computing W^, 
combined with different approximations for the functionals. 

Since the only purpose of this paragraph is to investigate the wavefunction 
part of the multideterminantal DFT, we only consider two very small systems, the 
He and the Be atoms, for which it has been possible to construct very accurate 
potentials v^ 1 |18| . In this way (i) we have essentially no approximation on the 
functional part and we can focus our attention on the effect of approximations on 
the calculation of ^f^, and (ii) we can produce accurate benchmark results to test 
our calculations. Of course, this is not the general procedure in which we are finally 
interested. The general procedure rather consists in using a given approximation 
(see Sec. I3.3|) for the functional E^. c [n] and the corresponding potential, and com- 
bining it with a wavefunction method to calculate to finally obtain the total 
ground state energy of the physical hamiltonian. This procedure has been followed 
with very promising results in the works reviewed in the next Sec. 0] 

Here, in order to illustrate the efficiency of approximate wavefunction meth- 
ods to treat we thus proceed as follows |18j . We first construct, for each 
fi, the model Hamiltonian using an accurate potential |18| and compute 
accurately its ground-state energy, E^ = (^^IH^I^^), at the multi-reference con- 
figuration interaction with singles and doubles (MRCISD) level. We then compute 
various approximate ground-state energies, Eg = | M | ^ ^) , by using approx- 
imate configuration inteaction (CI) type wave functions expanded into linear 
combinations of all the few Slater determinants generated from small orbital spaces 
S. The orbitals used are the natural orbitals of the Coulombic system calculated 
at the MRCISD level. The accuracy of the approximation for can be assessed 
by looking at the difference between Eg and E^ 

(3.17) AEg 1 = J5g - J5". 
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Figure 2. Ground-state energy differences AE£ = 
(^f^^H^l^^) where ^ is an accurate wave function and fyg are 
approximate wave functions generated from small orbital spaces 
S = Is, S = ls2s and S = Is2s2p, as a function of n for the He 
atom. 

The differences AEg are plotted as a function of /i in Fig. [21 for the He atom 
with the orbital spaces S = Is, S = ls2s and S = \s2s2p. One sees that, in the 
Coulombic limit, /x — > oo, the reduction of the orbital space leads to important 
errors in the energy. When /j, is decreased, i.e. when the interaction is reduced, the 
errors due to limited orbital spaces get smaller and smaller. For instance, at p = 1, 
using only the single-determinant wave function ^i sl leads to an error A£Jf s of less 
than 0.005 Hartree. 

The case of the Be atom with the orbital spaces S — ls2s and S — Is2s2p 
is reported in Fig. |3] Because of the near-degeneracy of the 2s and 2p levels, the 
inclusion of 2p configurations in the wave function is important, quite independently 
of the electron-electron interaction. Indeed, the difference E^ s2s — E^ s2s2p remains 
large for almost all /it's. On the contrary, the error of the calculation where the 2p 
orbitals are included, AE^ s2s2p , quickly falls off when p is decreased. Again, for 
fx = 1 for instance, the error AE^ s2s2p given by the few-determinant Cl-type wave 
function ^i s2s 2p 1S ^ ess than 0.005 Hartree. 

Therefore, the modification of the interaction enables to increase the accuracy 
of CI- type wave function calculations, or equivalently for a fixed target accuracy, 
decrease the effort of the calculation by reducing the orbital space. The crucial point 
for this effect to appear seems to be the reduction of the electron-electron interaction 
compared to the Coulomb interaction rather than the long-range character of the 
modified interaction. 

4. Some recent results: a short overview 

The idea of using the partially interacting model system of Eq. I|2.1|l has been 
explored with different techniques to compute the wavefunction and using 
different approximations for the functionals of Sec. 13.11 We review very briefly 
some of the corresponding results. 
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Figure 3. Ground-state energy differences = 
- ^\H^\^) where is an accurate wave 

function and are approximate wave functions generated from 
small orbital spaces S = ls2s and S — \s2s2p, as a function of fi 
for the Be atom. 



The LDA functional for E xc [n] has been combined with configuration interac- 
tion (CI) to handle the wavefunction ^ in Refs. [331 134 j ; the method has been 
applied to atoms and small molecules. Again using LDA for the short-range func- 
tional, in Ref. |35j it has been shown that long-range van der Waals forces in rare 
gas dimers can be well described by using second-order perturbation theory for ty^. 

A short-range GGA functional has been combined with the coupled-cluster 
(CC) method to describe "I^, with very good results for small molecules, both 
for the closed and the open shell cases |26l 127) : in particular, the results from 
multidctcrminantal DFT for small and medium basis-set sizes are better of both 
the pure DFT result and the pure CC result. In Ref. 36 , a different GGA func- 
tional |25l 118) has been used, and has been determined by multi-configuration 
self consitent field (MCSCF). The corresponding application to systems in which 
near-degeneracy effects play a major role is promising, although self-interaction 
errors in the functionals are still a problem in some cases |36| . 

Realistic / M (7*12) from Eqs. i|3.16|l have been generated for two-electron atoms 
in Ref. |29| : with a very simple approximation for w^ s (r 12) (inspired to the one 
used for the electron gas at /1 = 00), extremely accurate .E^cM have been obtained. 
This new strategy of producing functionals is still at a very early stage, and several 
steps are needed for its development. In particular, we need better approximations 
for Wefi(r 12), and the implementation of an efficient algorithm to couple the radial 
equations l|3.16|l with the method used to solve the model hamiltonian H^. 

In conclusions, changing the model Hamiltonian in DFT calculations seems 
to be a promising alternative to the construction of better exchange-correlation 
functionals for standard Kohn-Sham DFT. Of course, the method has been under 
development only in the last decade, and it still needs further improvement and 
investigation in many aspects. 
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